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Abstract. We present limiting shape results for a non-abelian variant of the abelian 
sandpile growth model (ASGM), some of which have no analog in the ASGM. One of 
our limiting shapes is an octagon. In our model, mass spreads from the origin by the 
toppling rule in Zhang's sandpile model. Previously, several limiting shape results 
have been obtained for the ASGM using abelianness and monotonicity as main tools. 
As both properties fail for our model, we use a new proof technique: in our main 
proof, we introduce several cellular automata to mimic our growth model. 
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1. Introduction 

We consider the following setup. Start with a pile of mass n > 1 at the origin of the 
rectangular grid 1/ , and mass /i < 1 at every other site, where n and h are real numbers. 
Mass may be moved by 'splitting piles', that is, one may take all the mass from one site, 
and divide it evenly among its neighbors. One can only split 'unstable' piles of mass at 
least 1, and one can not stop before all the piles have mass less than 1. We call T the - 
possibly random - set of sites where at least one split was performed. 

Will the mass spread over all of Z'^, or will T be a finite subset of lAl In the first case, 
how does it spread? In the last case, what is the size and shape of this set, depending on 
h and nl How do these answers depend on the order of splitting? 

This splitting game is related to abelian sandpile growth models. In fact, the 'splitting 
pile' rule is the same as the toppling rule in Zhang's sandpile model [10], a sandpile model 
that has not been considered as a growth model before. A notorious difficulty of moving 
mass by the Zhang toppling rule is the non-abelianness, that is, the end result depends 
on the order of topplings O |S] . 

In other studies of sandpile growth models [2||4l|8], ample use was made of the freedom, 
by abelianness, to choose some convenient toppling order. In this way, information could 
be derived about limiting shapes, growth rates and about whether an explosion occurs 
or not. The term 'explosion' is introduced in if an explosion occurs, then the mass 
from the origin spreads over all of lA, and every site topples infinitely often. For 'robust 
backgrounds', that is, values of h such that an explosion never occurs, one can examine 
the growth rate and existence of a limiting shape as ri — > 00. Roughly speaking, the 
growth rate is the radius of T as a function of n, and if the set T, properly scaled, 
converges to a deterministic shape as 71 —> 00, then that is the limiting shape. 



In [8], the main topic is the rotor router model. In the rotor router model, the mass 
consists of discrete grains, so that only discrete values of h are possible. All sites except 
the origin start empty, and at all sites there is a router which points to a neighbor. 
Instead of splitting, a pile that has at least one grain may give a grain to the neighbor 
indicated by the router. The router then rotates to point to the next neighbor, for a 
cyclic ordering of all the neighbors. It was proved that the limiting shape is a sphere, 
and the growth rate is proportional to n^^'^. The proof makes use of properties of Green's 
function for simple random walk. In this paper, we will demonstrate that the method of 
[5] can be adapted for our model, resulting in a growth rate proportional to n^/'* for all 



h < 0: see Theorem 3.4 



In the abelian sandpile model (for background, see [Il[9|), the mass also consists of 
discrete sand grains. Instead of splitting, a pile that consists of at least 2c? grains may 
topple^ that is, give one grain to each of its neighbors. Thus, 2d consecutive rotor router 
moves of one site equal one abelian sandpile toppling of that site. In 4J, making use of 
this equivalence, it was shown that for the abelian sandpile growth model with h — —oo, 
the limiting shape is a sphere. Moreover, it was proved that for h — 2d — 2, the limiting 
shape is a cube. In [2], it was proved that for all h < 2d— 2, the growth rate is proportional 
to n^/'^. 

All these proofs heavily rely on the abelianness, which is an almost routine technique 
for the abelian sandpile model, but fails to hold in our case. For instance, consider the 
case d = 1, n = A, and h = 0. If we choose the parallel updating rule that is common in 
cellular automata, in each time step splitting every unstable site, then we end up with 

. . . , 0, 1/2, 3/4, 3/4, 0, 3/4, 3/4, 1/2, 0, . . . 

However, if we for example choose to split in each time step only the leftmost unstable 
site, then we end up with 

...,0,1/2,1/2,7/8,3/4,0,3/4,5/8,0,... 

For arbitrary splitting order, it may even depend solely on the order if there occurs an 
explosion or not; see the examples in [3], Section 4.1. In this paper, we focus primarily 
on the parallel splitting order, but several of our results are valid for arbitrary splitting 
order. 

Another complicating property of our model is that unlike the abelian sandpile model, 
we have no monotonicity in h nor in n. In the abelian sandpile growth model, it is almost 
trivially true that for fixed h, T is nondecreasing in n, and for fixed n, T is nondecreasing 
in h. For our growth model however, this is false, even if we fix the parallel splitting order. 
Consider the following examples for c? = 1: For the first example, fix /i = 23/64. Then 
if n = 165/32 « 5.16..., we find that T is the interval [— 5, — 4, . . . , 4, 5]. However, 
if n = 167/32 « 5.22..., then T = [— 4, — 3, . . . , 3, 4]. For the second example, fix 
n = 343/64. Then if 21/64 « 0.33..., we find that T = [-5, ~4, . . . , 4, 5], but if 
h = 23/64 w 0.36 . . ., then T = [-4, -3, . . . , 3, 4]. 

One should take care when performing numerical simulations for this model. Since in 
each splitting, a real number is divided by 2d, one quickly encounters rounding errors 
due to limited machine accuracy. All simulations presented in this paper were done by 
performing exact calculations in binary format. 

This article is organized as follows: after giving definitions in Section [2] we present our 
results and some short proofs in Section [3j Our main result is that for h explosive and 
with the parallel splitting order, the splitting model exhibits several different limiting 
shapes as i — > oo. We find a square, a diamond and an octagon. These results are stated 
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in our main Theorem |3.H which is proved in Section |4] Sections [5] and [6] contain the 
remaining proofs of our other resuhs. Finally, in Section [7] we comment on some open 
problems for this model. 

2. Definitions 

In this section, we formally define the splitting model. While we focus primarily on 
the parallel order of splitting, some of our results are also valid for more general splitting 
order. Therefore, we give a general definition of the splitting model, with the splitting 
automaton (parallel splitting order) as a special case. 

For n G [0,oo) and h S (— (X),l), 77^ is the configuration given by 

,1 , , J n if X = 0, 
VnW - I if xgZ'^\0. 

For every ^ = 0,1,2,..., and for fixed n and h, rjt is the configuration at time t, and the 
initial configuration is rjo — rjl^. We interpret 77t(x) as the mass at site x at time t. 

We now describe how rjt+i is obtained from r]t, for every t. Denote by lAt — {x. : 
^(x) > 1} the set of all unstable sites at time t. St+i is a (possibly random) subset of 
Ut- We say that St+i is the set of sites that split at time t + 1. Then the configuration 
at time t + 1 is for all x defined by 

77t+i(x) = 7;t(x) (1 - lxe5t+i) + ^ ^ f7t(y)ly~x- 

ye5t+i 

The splitting order of the model determines how we choose St+i, given lAt, for every 
t. For example, if we have the parallel splitting order then we choose St+i = Ut for 
every t. In this case, we call our model the splitting automaton. Some of our results are 
also valid for other splitting orders. In this paper we only consider splitting orders with 
the following properties: At every time t, St+i is non-empty unless lAt is empty, and for 
every x that is unstable at time t, there exists a finite time t^ such that x € St+to- For 
example, we allow the random splitting order where St+i contains a single element of 
Ut, chosen uniformly at random from all elements of Ut- With this splitting order, at 
each time step a single site splits, randomly chosen from all unstable sites. This splitting 
order is valid because Ut increases slowly enough. Since only the neighbors of sites that 
split can become unstable, we have for any splitting order that Ut+i ^Ut + dUt, where 
with dX for a set X C Z"^, we denote the set of sites that are not in X, but have at least 
one neighbor in X. But when every time step only a single site splits, at most 2c? sites 
can become unstable, so that \Ut+i\ < \Ut \ + 2d. 

We are interested in the properties of 

7^= U 

0<t'<t 

all the sites that split at least once until time t, as well as 

T = lim Tt- 

For a fixed splitting order, we say that 77^ stabilizes if in the limit t — )■ 00, for every x, 
the total number of times that site x splits is finite. Note that if 77^ does not stabilize, 
then every site splits infinitely often. We also remark that in order to show that 77^ does 
not stabilize, it suffices to show that T = Z''. Namely, if T = Z'', then every site splits 
infinitely often. Otherwise, there is a site x and a time t such that x does not split at 
any time t' > t, but each of its neighbors y^, i = 1, 2, . . . , 2d, splits at some time ti > t. 
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But then at time maxi ti, x is not stable, because it received at least mass from each 
of its neighbors. Therefore, x must split again. 

We call stabilizable if 77^ stabilizes almost surely (The "almost surely" refers to 
randomness in the splitting order). As defined in 

Definition 2.1. The background h is said to be robust j/?7^ is stabilizable for all finite 
n; it is said to be explosive if there is a N'^ < 00 such that for all n > iV'\ ry^ is not 
stabilizable. 

Remark We expect the splitting model for every h to be either robust or explosive. 



independent of the splitting order (see Section 7.1 1. However, even for a fixed splitting 
order we cannot a priori exclude intermediate cases where the background is neither 
robust nor explosive. For example, since the splitting model is not monotone in n, it 
might occur for some h that for every n, there exist ni > n2 > n such that ry^^ is 
stabilizable, but 77^^ is not. 

Finally, we give our definition of a limiting shape. In ^IHIl], limiting shape results 
were obtained for T in the limit n — >■ 00, for robust background. In this paper however, 
we present limiting shape results for Tt with n fixed, parallel splitting order and explosive 
background. We study the limiting behavior in t rather than in n. Accordingly, we have 
a different definition for the limiting shape. 

Let C denote the cube of radius 1/2 centered at the origin {x e M'' : max^ Xi < 1/2}. 
Then x + C is the same cube centered at x, and by V + C we denote the volume 
Uxev(^ + C). 

Definition 2.2. Let Vt, with t = 0,1,... be a sequence of sets in 1'^ . Let S be a 
deterministic shape in M"*, scaled such that maxa,.{x S S} = 1. Let and S"^ denote 
respectively the inner and outer e -neighborhoods ofS. We say that S is the limiting shape 
ofVt, if there is a scaling function f(t), and for all e > there is a such that for all 
t > t\ 

SeC/(<) (Vt + C)CS^ 
// S is the limiting shape ofTt, then we say that S is the limiting shape of the splitting 
automaton. 



3. Main results 



We have observed - see Figure [T]- that varying h has a striking effect on the dynamics 
of the splitting automaton. For large values of h, Tt appears to grow in time with linear 
speed, resembling a polygon, but which polygon depends on the value of h. Figures [T] and 
|3] support the conjecture that for the parallel splitting order, as t ^ 00, there are many 



possible different limiting shapes depending on h. Our main result, Theorem 3.1 is that 
there are at least three different polygonal limiting shapes, for three different intervals 
of h. For small values of /i, the splitting model stabilizes; see Theorem |3.2[ In between, 
there is a third regime that we were not able to characterize. It appears that for h in 
this regime, Tt keeps increasing in time, but does not have a polygonal limiting shape. 
We comment on this in Section [71 
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Figure 1. The splitting automaton for different values of h. "Warmer" 
color is larger mass; dark yellow, orange or red cells have mass > 1. Top 
row: h = 47/64 « 0.734; h = 1495/2048 « 0.73 and h = 727/1024 w 
0.71, each with n — 8 and after 50 time steps. Bottom row: ft, = 1/2 
and n = 256; h = 511/1024 w 0.499 and n = 2048; h = and n ^ 2048, 
each after the model stabilized. 




Figure 2. The diamond D, the square Q, and the octagon O. 
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Let D be the diamond in M'^ with radius 1 centered at the origin; let Q be the square 
with radius 1 in centered at the origin, and let O be the octagon in R-^ with vertices 
(0,1) and (§,§), and the other six vertices follow from symmetry. See Figure [2]for these 
shapes. 

Theorem 3.1. 

(1) The limiting shape of the splitting automaton on , for 1 ~ < h < 1 and 
n > 1, is the diamond D. The scaling function is f{t) = 

(2) The limiting shape of splitting automaton on 1? , for 7/10 < h < 40/57 and 
A — Ah < n < 16 — 20h, is the square Q. The scaling function is f(t) — |. 

(3) The limiting shape of the splitting automaton on 1? , for hjl < h < 13/18 and 
n — 3, is the octagon O. The scaling function is f{t) = 

In the abelian sandpile growth model, h = 2d — 1 is the only possible explosive back- 
ground value, and our proof of part 1 also works for that situation. However, the second 
two parts have no analog in the abelian sandpile growth model. 

Our proof uses the method of mimicking the splitting automaton with a finite state 
space cellular automaton; we consider our three explicit results as introductory examples 
for this method. We expect that with this method, many more limiting shape results can 
be obtained. 

Next, we characterize several regimes of h for the splitting model on Z"^. 

Theorem 3.2. In the splitting model on , 

(1) The background is explosive if h > 1 — 

(2) The background is robust if h < ^, 

(3) In the splitting automaton, for d > 2, there exist constants < 1 — ^^"^^ such 
that the background is explosive if h > Cd- 

We give the proof of the first part here, because it is a very short argument. The proof 
of parts 2 and 3 will be given in Section [5] where we give the precise form of Cd- We 
do not believe that this bound is sharp. From the simulations for d — 2 for example, 
a transition between an explosive and robust regime appears to take place at h = 2/3, 
while in our proof of part 3, C2 ~ 13/19 w 0.684. 

Proof of Theorem \3.S\ part 1 

If n > 1, then the origin splits, so then T is not empty. Now suppose that T is a finite 
set. Then there exist sites outside T that have a neighbor in T- Such a site received at 
least ^, but did not split. For ft- > 1 — 5^, this is a contradiction. □ 

Note that Theorem |3 . 2 1 does not exclude the possibility that there exists, for d fixed, a 
single critical value of h that separates explosive and robust backgrounds, independent of 
the splitting order. We only know this in the case d — \^ for which the first two bounds 
are equal. 

We give another result that can be proved by a short argument: 

Theorem 3.3. In every splitting model on 1/ , for every n > 1 and h > 1 — ^; */ the 
model stabilizes then T is a d-dimensional rectangle. 

Proof. Because n > 1, T is not empty. Suppose that T is not a rectangle. Then, as is 
not hard to see, there must exist a site that did not split, but has at least two neighbors 
that split. Therefore, its final mass is at least ft- + g, but strictly less than 1. This can 
only be true for ft < 1 — i. □ 
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Figure 3. The splitting automaton after 50 time steps, with n = 3 and 
h=l= 0.75 (diamond), = § « 0.72 (octagon), and /i = |f| « 0.701 
(square) . 

In the case of the paraUel sphtting order, we additionally have symmetry. Thus, for 
/i>l— i,Tisa cube. We remark that the above proof also works for the abelian 
sandpile model, thus considerably simplifying the proof of Theorem 4.1 (first 2 parts) in 

El- 

Our final result gives information on the size and shape of T when h < 0. This theorem 
is similar to Theorem 4.1 of [8 . Let B^. denote the Euclidean ball in R'^ with radius r, 
and let ujd be the volume of Bi. 

Theorem 3.4. 

(1) (Inner bound) For all h < 1, 

with ci = (1 — h)^^/'^ , r — (-^)i/'^ and Ci a constant which depends only on d; 

(2) (Outer hound) When h <0, for every e > 0, 

with c'l — — e ~ h)^^/'^, r — i'^)^^'^ o,nd C2 a constant which depends only on 
e, h and d. 

4. Limiting shapes in the explosive regime 

In this section, we will prove Theorem |3.1[ Each part of this theorem is stated for h 
in a certain interval, and the first two parts for n in a certain interval. That means the 
theorem is stated for uncountably many possible initial configurations. However, we will 
show that we do not need to know all the exact masses to determine 7t for a certain t. 
For each part of the theorem, we will introduce a labeling of sites, using only a finite 
number of labels. It will suffice to know the labels of all sites at time t, to determine 7t' 
for all t' >t. 

We will see in each case that the time evolution in terms of the labels is a lot more 
enlightening than in terms of the full information contained in rjt- In each case, we can 
identify a certain recurrent pattern of the labels. Our limiting shape proofs will be by 
induction in t, making use of these recurrent patterns. 

The label of a site at time t will depend on its own label at time t—1 plus those of its 
neighbors at time t—1. We will specify the labels at t = 0, and the transition rules for 
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the labels. In other words, for each part of the theorem we will introduce a finite state 
space cellular automaton, that describes the splitting automaton for certain intervals of 
h and n in terms of the labels. 

A cellular automaton is defined by giving its state space S, its initial configuration 
^0, and its transition rules. By ^t(x), we denote the label of site x at time t. The state 
space will consist of a finite number of labels. The transition rules specify how the label 
of each cell changes as a function of its own current label and those of its neighbors. A 
cellular automaton evolves from the initial configuration in discrete time; each time step, 
all cells are updated in parallel according to the transition rules. 

We use the following notation for the transition rules. Let s' and s, si, S2, . . . denote 
labels in S (not necessarily all different). By 

S ® Si,S2, ■ ■ ■ ,S2d > s', 

we mean that if a cell has label s, and there is a permutation of the labels of its 2d 
neighbors equal to {si,S2, ■ ■ ■ , S2d}, then the label of this cell changes to s' . By *, we will 
denote an arbitrary label. For example, if we have a transition rule 

S © Si, Si, *, . . . , * !• s', 

then the label of a cell with label s will change into s' if at least two of its neighbors have 
label si, irrespective of the labels of the other neighbors. 

We first give the proof for part 1, which is the simplest case. In fact, in this case 
the splitting model is equivalent to (l,d) bootstrap percolation: Tt+i is the union of 7t 
with all sites that have at least one neighbor in Tt- The proof we give below will seem 
somewhat elaborate for such a simple case. That is because we use this case to illustrate 
our method of labels and cellular automata. 

We will need the following observation, which can be proved by induction in t: 

Lemma 4.1. We call x an odd site if^^Xi is odd, otherwise we call x an even site. 
Then in the splitting automaton, even sites only split at even times, and odd sites only 
at odd times. 

Proof of Theorem \3.1\ part 1 

We begin by defining the diamond cellular automaton. 
Definition 4.2. 

The diamond cellular automaton has state space {e,h,u\. We additionally use the 
symbol s to denote e or h. In the initial configuration, every cell has label h, only the 
origin has label u. The transition rules are: 

(1) h® s,...,s — >h 

(2) s ® u, . . . ,u — > u 

(3) M © s, . . . , s — > e 

(4) /i© u, * — > u 

The above set of transition rules is sufficient to define the diamond cellular automaton, 
because, as we will demonstrate below, other combinations of cell and neighborhood labels 
do not occur. 

Let Qt, the growth cluster of the cellular automaton, be the set of all cells that do 
not have label h at time t. We will first prove the limiting shape result for Qt, and then 
demonstrate that ifl— ^</i<l, then Qt is the same set as 7t U dTt for every t. 

Let T>r be the diamond {x g Z"^ : \xi\ < r}. To prove the limiting shape result, 
we will show by induction in t that Qt —'Dti so that the limiting shape of CJ( is D, with 
scaling function f{t) = i. 
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Our induction hypothesis is that at time t, S^t is as follows (see figure [4|: all sites 
X € Vt have label u if Xi — t) mod 2 = 0, and label e otherwise. All other sites have 
label h. If this claim is true for all t, then we have Gt = 2^t- 
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Figure 4. The induction hypothesis for the diamond cellular automaton 
aX t = 2. Labels not shown are h. 

As a starting point, we take t — 0. At that time, the origin has label u, and all other 
cells have label h. Therefore, the hypothesis is true at t = 0. 

Now suppose the hypothesis is true at time t. Then all sites with label e have 2d 
neighbors with label m, therefore by the second transition rule they will have label u at 
time t + 1. All sites with label u have 2d neighbors with label e or ft, therefore by the 
third rule they will have label e at time t + 1. All sites in Vt+i \ Vt have label ft and at 
least one neighbor with label u. Therefore, by the fourth rule they will have label u at 
t + 1. Other labels do not change, by the first rule. This concludes the induction, and 
moreover shows that our set of transition rules is sufficient to define the diamond cellular 
automaton. 

Finally, we show that ifl — ^ < h < 1, then Qt is the same set as 7t U dTt- To compare 
the configurations ryt and ^t, we give a mapping 

Md ■ {e,ft, M, s} -^I, 

where I is the set of intervals {[a, b) : a < b, a,b E [0, oo]}, that maps the state space of 
the diamond cellular automaton to the mass values of the splitting automaton: 

Md{e) = 

Md{h) = ft 

Md{s) = [0,1) 

Md{u) = [l,oo) 



For a fixed diamond cellular automaton configuration ^, define 
MI = {7j: r/(x) e Xd(e(x)), for all x e Z'*}. 

With this mapping, the initial configuration of the splitting automaton 77^ is in A^^", 
for all n > 1 and 1 — ^ < ft < 1. We will show by induction in t that % is in M.^ for 
all t. Suppose that at time t, rjt is in A^^'. 

We check one by one the transition rules: 

• (rule 1) If in the splitting automaton a site has mass ft, and none of its neighbors 
split, then its mass does not change. This is true for all ft < 1. 

• (rule 2) If in the splitting automaton a stable site has 2d unstable neighbors, then 
its receives at least 2d times therefore its mass will become at least 1. This 
is true for all ft. 
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• (rule 3) If in the splitting automaton a site has mass at least 1, then it splits. If 
none of its neighbors splits, then its mass will become 0. This is true for all h. 

• (rule 4) If in the splitting automaton a site with mass h has at least one neighbor 
that splits, it receives at least Therefore, it will become unstable only if 

l-h<h<l. 

Therefore, if ?jt is in M^^\ 1 — ^ < ^ < 1 and n > 1, then ijt+i is in TM^'^^ This 
completes the induction. 

Finally, by the following observations: 

• only the label u maps to mass 1 or larger, so a site is in 7t if and only if it has 
had label u at least once before t, 

• the label of a site changes into another label if and only if at least one neighbor 
has label u, 

• if a site does not have label h at time t, then it cannot get label h at any time 
t' > t, 

we can conclude that Qt of the diamond cellular automaton is the same set as 7t U dTt 
for the splitting automaton with 1 — ^ < h < 1 and n > 1. □ 

We now give the proofs of the remaining two parts; note that in these next two proofs, 
we are in dimension 2. We will need more elaborate cellular automata, in which there 
are several labels for unstable sites. For example, it is important to know whether the 
mass of a splitting site is below or above 4(1 — h): if its neighbor has mass h then in the 
first case it might not become unstable, but in the second case, it will. 

Proof of Theorem \3.1\ part 2 

We begin by defining the square cellular automata. 
Definition 4.3. 

The square cellular automaton has state space {e,h,p,m,m\c,d}. We additionally 
use the symbol s to denote a label that is e, h or p. In the initial configuration, every cell 
has label h, only the origin has label d. The transition rules are: 
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3 d, d, c,p — > 


to' 



The above set of transition rules is sufficient to define the cellular automaton, because, 
as we will demonstrate below, other combinations of cell and neighborhood labels do not 
occur. 

Recall that the growth cluster Qt is the set of all cells that do not have label h at 
time t. To prove the limiting shape result for the growth cluster of the square cellular 
automaton, we use induction. Let Cr G be the square {{i,j) ■ \i\ < r, \j\ < r}. Let 
be the following configuration (see Figure [s]): 

• All sites in C^-i have label c. 

• The labels in Cr \ Cr^i are d, if [i — j) mod 2 = 0, and e otherwise. 

• The labels outside Cr are p if they have a neighbor with label e, and h otherwise. 
Our induction hypothesis is that for every even t, = The initial configuration 

of the square cellular automaton is Co- Now suppose that at some even time t, — Ct/2- 
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Figure 5. The configuration Cr, used in the induction hypothesis for 
the square cellular automaton, for r ~ 2. Labels not shown are h. 



Then by using the transition rules, one can check that at time S,t+2 will be Ci/2+1 ~ 
C(t+2)/2- This completes the induction. 

The shape result now follows: for every i, Ct/2 Q Gt Q Ct/2+1, so that the limiting 
shape of Gt is the square Q, with scaling function f{t) — |. 

Finally, to show that Gt for the square cellular automaton is equal to 7t U dTt for 
the splitting automaton with 4 ~ Ah < n < 16 — 20h and 7/10 < h < 40/57, we give a 
mapping 

A^s : {e, ™, ITT-', c, d, s} — ?► I, 

that maps the state space of the square cellular automaton to the mass values of the 
splitting automaton: 



Ms{e) 


= 


MM 


= h 


Ms{s) 


= [0,1) 




= [h+\,l) 


Ms{m) 


= [1,4-4/1) 


Ms{rn') 


= [0,12-15/i) 


Ms{c) 


= [0, 16 - 20/i) 


Ms{d) 


= [4 - 4/1, 16 - 20/1) 



For all h < 3/4, these intervals are nonempty, moreover, 4 — 4/i > 1. 

With this mapping, one may check that is in By induction in t, we will show 

that rjt is in for all t. Suppose at time t, rjt is in A^|*. 

We check one by one the transition rules: 

• (rule 1) If in the splitting automaton a site has mass h, and none of its neighbors 
split, then its mass does not change. This is true for all h < 1. 

• (rules 2-5) If an unstable site splits, then by Lemma 4.1 its neighbors do not split. 
Therefore, it will become empty. If a cell has split at least once, then from that 
time on it cannot receive sand from more than 4 neighbors before splitting itself. 
Therefore, no cell that split at least once can gain mass greater than 16 — 2Qh. 
This is true for all h < 1. 

• (rule 6) /i + i [4 - 4h, 16 - 20h) ^ [1, 4 - 4/i). This is true for all h < 1. 
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• (rule 7) /i + |[1,4 - Ah) [h + 1/4, 1). This is true for all h<l. 

• (rules) [/i+3,l) + i[l,4-4/i) + i[0,12-15/i) ^ [/i+f,6-^). We have that 
[h + f , 6 - ^) C [4 - 4/i, 16 - 20h) if 13/20 <h< 40/57. 

• (rule 9) h + ^[1,4 - Ah) -> [/i + ^,2 - h). We have that [/i + i,2 - /i) C 
[4 - Ah, 16 - 20h) if 7/10 < /i < 14/19. 

• (rule 10) h+ l[A- Ah, 16 - 20h) + l[l,A- Ah) [|,5 - 5h). We have that 
[|, 5 - 5h) C [4 - Ah, 16 - 20/i) if 11/16 <h< 11/15. 

• (rule 11) i[4-4/i,16-20/i) + i[0, 16-20/i) ^ [2 - 2/i, 12 - 15/i). We have that 
[2 - 2h, 12 - 15h) C [0, 12 - I5h) if /i < 1. 

Therefore, if r/t is in Mf , A-Ah<n<m-20h and 7/10 <h< AO/57, then r^t+i is 
in Adi**^ . This completes the induction. 
Finally, by the following observations: 

• the labels m and d map to an interval in [1, oo), so a site is in Tt if it has had 
label m or c? at least once before t, 

• the label of a site with label h changes into another label if and only if at least 
one neighbor has label m or d, 

• if a site does not have label h at time t, then it cannot get label h at any time 
t' > t, 

we can conclude that ii A ~ Ah < n < 16 - 20h and 7/10 < h < 40/57, then Gt for the 
square cellular automaton is the same set as Tt U dTt- □ 

For the final part, we first perform 8 time steps in the splitting automaton before we 
describe its further evolution as a finite state space cellular automaton. Otherwise, we 
would need many more labels and transition rules. 

Proof of Theorem \3.1\ part 3 

We begin by defining the octagon cellular automaton. 
Definition 4.4. 

The octagon cellular automaton has state space {e, h,p, m, d, d', d\, c, c' , q, q'}. We ad- 
ditionally use the symbol s to denote a label that is e, h or p, and the symbol u to denote 
a label that is any of the other. 

The initial configuration is given in the table below. We show only the first quadrant 
(left bottom cell is the origin). The rest follows by symmetry. All labels not shown are h. 
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P 






e 


c 


e 


c 


e 


d! 






c 


e 


c 


e 


c' 


e 


P 





The transition rules are: 



(1) hQ 


s,s,s,s — 


-> h 


(6) hQ 


3 d\, s,s,s — 


-> TO 


(2) uQ 


— 


e 


(7) hQ 


^q,s,s,s — 


> d 


(3) hQ 


m, s, s, s - 


->P 


(8) hQ 


3 <7, TO, s, s — 


d\ 


(4) ha 


^ d,s,s,s — 


m 


(9) hQ 


3 q' , m,s,s ' 


^ d\ 


(5) hQ 


5 d' , s,s,s — 


— > TO 


(10) hQ 


3 TO, d', s, s ~ 


^d' 
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(11) h®d',d',s,s — > d' (17) eS)q,c,dl,s — >c 

(12) p © d\, 171,0,3 — > q' (18) e ® q', c, d\,s — > c 

(13) p(Bm,m,c,s — >q (19) e®u,u,u,u — >c 

(14) p (B d,m,c, s — > q except in the following case: 

(15) e®q',c,d',s — > c (20) e®m,q,q,c — > c' 

(16) e®dl,d\,c',s — > c 

Our proof that the growth cluster of the octagon cellular automaton has O as limiting 
shape, is by induction. We will show that there is a pattern that repeats every 10 time 
steps. To describe this pattern, we introduce two subconfigurations that we call 'tile' and 
'cornerstone'. They are given in Figure |6] We say the subconfiguration is at position 
{x,y), if the left bottom cell has coordinates (x,y). 
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Figure 6. A 'tile' (left), and a 'cornerstone' (right). 

To specify a configuration, we will only specify the cells x = {x, y) with y > x > 0; 
the rest follows by symmetry. 

Definition 4.5. We define the configuration x*, with i G {0, 1, . . as follows: 

• there is a cornerstone at position {b{i + 1), b{i + 1)), 

• for every j = 0, . . . ,i, there is a tile at position (5(i — j),7 + 5i + j), 

• the leftmost tile is different, namely, xH^i 7 + 6i) = c' , 

• for every cell {x, y) such that it is not in a tile or cornerstone, hut every directed 
path from {x,y) to (0,0) intersects a tile or cornerstone, x^{x,y) — h, 

• for every other cell (x, y), X^{^^ 2/) = c if {x + y) mod 2 = 0, and c otherwise. 

Our induction hypothesis now is: In the octagon cellular automaton, at time 5 + 10«, 
Cs+ioi — X*- III words this says that every 10 time steps, an extra tile is formed. 

The hypothesis can be verified for i < 6, by performing 65 time steps starting from the 
initial configuration ^o- We show the results of this computation in Figure [TJ generated 
by a computer program of the octagon cellular automaton. 

Suppose now that the hypothesis is true at time 5 + lOi, with z > 6. We will construct 
^5_i_io(i+i) by performing 10 time steps starting from x*- Observe that the label of a cell 
at time t + 10 depends only on its own label and that of all cells in x + I?io; we call this 
set of cells the '10-neighborhood' of x. By the definition of x*, we have that for every 
i > b and for every x, there exists a cell y such that the labeling of the 10-neighborhood 
of X in X* is identical to that of the 10-neighborhood of y in x^- Therefore, the label of 
X in x*^^ will be identical to that of y in x^- Thus, we can construct ^5+10(^+1) from x\ 
and we find that indeed, if £,5+wi = X* then ^5+10(^+1) = x'^^- 

Since Qt is nondecreasing in t, we have that for every t there is an i such that Q^+iai C 
Qt ^ ^^5+10(1+1)- The radius of Q^+wi is 9-|-6i. This means that every 10 time steps, the 
radius increases by 6. We conclude that the limiting shape of Qt is the octagon O, with 
scaling function f{t) = ^. 
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Figure 7. The octagon cellular automaton after 5, 15, 25, 55, 65 and 
135 time steps, or equivalently, the configurations x", x^, X^j X^; X^ and 
X^^- Black = e, yellow = c, dark yellow = c', light blue — p, dark blue 
= fi,, light green = m, orange = d', red = d\, dark red = g'. 
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Finally, we prove that Gt of the octagon cellular automaton is equal to Tt+s U dTt+s, if 
n = 3 and e [5/7,13/18). 
We give a mapping 

Mo ■ {e, m, d, d' , d\, c, c', q, q' , s, u\ — > I, 

that maps the state space of the octagon cellular automaton to the mass values of the 
splitting automaton: (if 1/2 < h < 946/1301 « 0.727 then all intervals are nonempty, 
and moreover, the labels e, h and p map to an interval in [0, 1), while all other labels 
map to an interval in [1, cx))): 



Mo{e) 


= 


Mo{h) 


= h 


Moiv) 




Mo{s) 


= [0,1) 


Mo{m) 


= [1,4-4/1) 


Moid) 


= [4 - Ah, 16 - 20/i) 


Moid') 


= [^ + ^,16-20/^) 


Moidl) 


= [1 + ^,16-20/.) 


Moiq) 


= [1 + /i, 60 - 80/i) 


Moiq') 




= [^ + \m~m 


Moic) 


= [l,60- 80/i) 


Moic') 


= [1 + ^,60-80/i) 


Moiu) 


= [l,60- 80/i) 



For every x, if n = 3 and /i e [5/7,13/18) then r/g e A^|«. This can be verified by 
tedious, but straightforward inspection: In Figure [8] we give the configuration at t = 8 
for the splitting automaton with n — i and h E [5/7, 13/18). 
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Figure 8. rjg, multiplied by 4^ — 65536, for the splitting automaton 
with n — 3 and h £ (5/7, 13/18). Masses not shown are 65536/i. 

Next, we will prove by induction in t that rjt+s is in M^' for all t. We assume that for 
some t, ?7t+8 is in . By examining every transition rule, we can then show that as a 
consequence, r]t+g is in A^|'+\ 

We check one by one the transition rules: 
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(rule 1) If in the splitting automaton a site has mass h, and none of its neighbors 
split, then its mass does not change. This is true for all h < I. 



(rule 2) If an unstable site splits, then by Lemma 4.1 its neighbors do not split. 
Therefore, it will become empty. 

(rule 3) /i + |[1,4 - Ah) [h + 1/4, 1). This is true for ah h < I. 
(rule 4) /i + i [4 - Ah, 16 - 20h) ^ [1, 4 - 4/i). This is true for edl h < 1. 
(rule 5) /i+ j[| + f^,16-20/i) ^ + 4 - 4/i). We have that [^ + ^^,4- 
Ah) C [1, 4 - 4/i) if 1 > /i > 29/42 06905. 

(rule6) h+m+'f,16-20h) -> [| + ^,4-4/i). We have that [i + ^,4-4/i) C 
[1,4- Ah) if 1 > /i > 2/3. 

(rule?) h+\[l+h,60-8Qh) [^ + l,15-19h). We have that [f^ + i, 15-19/i) C 
[4 - Ah, 16 - 20h) if 0.7143 w 5/7 < /i< 1. 

(rule 8) h+\[l + h,60 - 80h) + \[1,A-Ah) ^ + i,16- 20/i). This is true 
for aU h <1. 

(rule 9) /i+i[^ + |,60- 80/i) + |[l,4-4/i)^ [f^ + ^, 16 - 20/i). We have 
that [f^ + ^, 16 - 20h) C + i, 16 - 20h) iil>h> 5/2. 
(rule 10) h + + §,16 - 20h) + ^[l,A-Ah) + 5h). We have 

'j^+'i,5-5h)C[f 



that + ^,5-5h) + |,16- 20/i) if 1/2 < /i < 11/15 w 0.7333 



8 ^ 16 ' 



• (rule 11) /i+i[f^ + |,16-20/i) ^ + ^, 8- 9/i). We have that 
9/i) C + 1^ 16 - 20h) if 1/2 < /i < 8/11 « 0.7273. 

• (rule 12) [h + 1,1) + i[f + 5,16 - 20h) + ^[1,4 - Ah) + i[l,60 - 80/i) ^ 

+ |, 21 - 26/i). We have that + i, 21 - 26/i) C + 1, 60 - 80h) if 
/i < 13/18 « 0.7222. 

• (rule 13) [h+l,l) + ^[l,A-Ah) + i[l,60- 80/i) ^ [1 + h,18-22h). We have 
that [1 + h,l8- 22h) C[l + h,60- 80h) if h < 21/29 « 0.7241. 

• (rule 14) [/i+i,l) + i[l,4-4/i)i[4-4/i, 16-20/i) + i[l,60-80/i) ^ [|,21-26/i). 
We have that [|, 21 - 26/i) C [1 + /i, 60 - 80/i) if /i < 13/18 w 0.7222. 

• (rule 15) i[^ + |,60 - 80/i) + i[l,60- 80/t) + i[f + |, 16-20/1) ^ [fj^ + 
^,34 - A5h). We have that [^^f + ^,34 - A5h) C [1,60 - 80h) if 0.6829 « 
28/41 <h< 26/35 w 0.7429. 

• (rule 16) i[f^ + i,16-20/i) + i[l + |,60-80/i) ^ + i, 23 - 30/i). We have 
that + i 23 - 30/i) C [1, 60 - 80h) if 2/3 <h< 37/50 w 0.7400. 

• (rule 17) i[l+/i,60-80/i) + i[l,60-80/i) + i[5^ + i,16-20/i) ^ [f| + |, 34-45/i). 
We have that [f| + |, 34 - A5h) C [1, 60 - 80h) if 2/3 < /i < 26/35 w 0.7429. 

. (rule 18) + 1,60 - 80/1) + i[l,60- 80/i) + i[f + i, 16- 20/^) ^ + 

i|,34 - A5h). We have that + if, 34 - 45/i) C [1,60 - 80h) if 0.6341 « 
26/41 < /i < 26/35 « 0.7429. 

• (rule 19) If in the splitting automaton a cell is empty, and all its neighbors split, 
then it gets mass at least 1. Since no cell has mass exceeding 60 — 8Qh, that is 
the maximum mass that an empty cell can get. 

• (rule 20) |[1,4- 4/i) + ^[1 + /i, 60 - 80/i) + i[l,60 - 80/i) ^ [1 + |,46- 61/i). 
We have that 46 - 6lh < 60 - 80h if /i < 14/19 « 0.7368. 

In summary, all the rules are valid if h £ [5/7, 13/18]. Therefore, if rjt+g is in A^|* and 
5/7 < h < 13/18, then rjt+g is in This completes the induction. 

Finally, by the following observations: 



• all sites in Qq are in Ts U dTs, 
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• only labels denoted as u map to values in [1, oo), so for all i > 8, a site is in 7t if 
and only if it is in Ts or it has had a label denoted as u at least once before t, 

• the label of a site with label h changes into another label if and only if at least 
one neighbor has a label denoted as w, 

• if a site does not have label h at time t, then it cannot get label h at any time 
t' > t, 

we can conclude that if n = 3 and h E [5/7,13/18], then Qt for the octagon cellular 
automaton is the same set as Tt+a U dTt+a- 

□ 



5. Explosive and robust regimes 
In this section, we prove parts 2 and 3 of Theorem |3. 2 1 



Proof of Theorem \3.S\ part 2 

We will prove that for all n, all /i < 1/2 and alH, |7t| < jj^-h^ where by |T| we denote 
the cardinality of a set T C Z*^. It follows that 

(5.1) \T\ < 



so that for all h < 1/2, the background is robust. 

Let mo be the total mass in 7t U dTt at time 0, and let rrit the total mass in 7t U dTt 
at time t. We have 

mo=n + h\Tt\+h\dTt\. 

At time t, TtU dTt contains a total mass of at least ^ times the number of internal 
edges in 7t U 97*. Namely, consider a pair of sites x and y connected by an internal edge. 
Each time that one of them splits, a mass of at least ^ travels to the other one. 

The number of internal edges in 7* U dTt is at least d\T\- We demonstrate this by the 
following argument: Fix an ordering for the 2d edges connecting a site to its 2c? neighbors, 
such that the first d edges of the origin are in the same closed orthant. For every site x 
in Tt, all its edges are in 7t U dTt- If for every x g 7* we count only the first d edges, then 
we count each edge in 7t U dTt at most once, and we arrive at a total of d\Tt\- 

Therefore, at least a mass of d^\Tt\ = ||7t| remains in 7t U dTt- Moreover, since 
the sites in dTt did not split, the mass h at every site in dTt also remains in 7t U dTt- 
Therefore, we have 

mt > l\Tt\ + h\dTt\. 

Finally, we note that since up to time t no mass can have entered or left T U dTt i we 
have mo = mt- Putting everything together, we find 

n + h\Tt\ + h\dTt\ > ^\Tt\ + h\dTtl 
from which the result follows. □ 

Proof of Theorem \3.^ part 3 

We first give the proof for d > 3. 

First, we need some notation. Denote by Vr C Z"^' the diamond Vr — {x ■- Xi < r}, 
and by £r C the layer {x : J^i^i — '^l- Denote dk = {k,k,---k) G - Let 
rfe,o — dk, and for i — 1 . . . d, let ^ be the set of sites in C^k+i that have i nearest 
neighbors in T^^i-i- Observe that for every i, F^ ^ is not empty, and that Fj, ,^ — d^+i- 
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For example, in dimension 3: di = (1,1,1), Ti i = {(1, 1, 2), (1, 2, 1), (2, 1, 1)}, Ti 2 = 
{(1, 2, 2), (2, 1, 2), (2, 2, 1)},^ and T.^^ = ^2 = (2, 2, 2). 

Let = Eta ^rf = (4^' ^'^'^ ^* '^^^'^^^ qd+Pdh* = 2d{l - h*), 

so that 

(7d - 2(i 

n = 

Pd + 2d 

Finally, we define 

= max{l - 2'^*^' 
and remark that < 1 — ^^'^^ ; with equality only in the case d — 2. 
We will prove the following statement: 

Lemma 5.1. In the splitting automaton on Z'', if h > C'^, and n > 2d(l — h), then for 
every k >\, at time dk + 2, the sites in Fj, 1 split. 



Theorem [3^2] part 3 follows from Lemma |5 . 1 1 combined with Theorem |3.3| Lemma |5.1 
tells us that for every r, there is a site on the boundary of the cube {a; : max^ Xi < r} 



that splits at least once. By Theorem 3.3 and by symmetry, all sites in this cube split at 
least once. Therefore, limt_j.oo Tt = 1 ■ 

Proof. Note that at time t, no site outside Vi can have split yet, so if a site in Ffc ^ splits 
at time dfc + j + 1, it does so for the first time. We will show that in fact the sites in 
Ffc i do split at time d/c + i + 1. Since we take the parallel splitting order, by symmetry, 
all sites in F^^j distribute the same mass when they split; denote by m(fc, z) the mass 
distributed in the first split of a site in F^^^. 

We will prove the lemma by induction. For fc = 0, the lemma is true, because we chose 
n large enough. Now suppose it is true for some value fc. Then at time dfc + 2, the sites 
in Lfc.a receive ^m(fc, 1) > because they each have two neighbors in F^^i. Since they 
did not split before, their mass is now at least /i + 3 . Therefore, if > 1 — 3 , they split 
at time dfc + 3. This condition is fulfilled because > 1 — 

Continuing this reasoning, we find for all i = 2, . . . , d — 1 that at time dfc + i + 1, the 
sites in Ffe^j split, because each site in F^^i has i neighbors in Ffc_i_i, so it receives mass 
^'m{k,i — 1) > ^ > ^. We calculate, using that for i — 2,...,d, we have m(k,i) = 
h + i^m{k, i — 1), 

m(0, d) = qdm{0, 1) + pdh > qa. + pdh. 

Recall that F^^^ = Ffe+i.Q. Therefore, at time d{k + 1) + 1 the sites in F^+i i receive 

mass "^^2(f'' ■ If + "'2^'^"' — I' then the sites in F^+i^i split at time d{k + 1) + 2. This 
condition is fulfilled if /i > C^. This completes the induction. Therefore, in Theorem 3.2 



part 3, for d > 3 we can take Cd = C'd. □ 

For d — 2, C'^ is equal to 1 — — 0.7. Li this case, we can take C2 = 13/19 = 

0.684...: 

Proposition 5.2. In the splitting automaton with d = 2, the background is explosive if 
h > 13/19. 



Proof. In the proof of Theorem 3.2 part 3, we have proved that at time 2fc + 2, sites 
(fc + 1, fc) and (fc, fc + 1) split, making use of the fact that at time 2fc + 1, site (fc, fc) splits. 
We did not take into account that more sites in £k might split at time 2fc + 1. 

We now choose n > 64 — 84/i, so that at i = 3, sites (0,2), (1,1) and (2,0) split 
for the first time. We prove by induction that ii h > 13/19, then at time 2fc, the sites 
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(fc — 1, fc + 1) and {k + 1, k ~ I) have mass at least 1, and site (fc, k) has mass at least 
2 — h > 1. With our choice for n, this is true for /c = 1. Now suppose the hypothesis is 
true for k. This impUes that at time 2k + 1, the sites sites (fc, fc + 1) and (fc + 1, fc) have 
mass at least ^ + f- If this is at least equal to 4 — Ah, then we obtain the induction 
hypothesis for fc + 1. Solving ^ + | > 4 - 4/i gives h > 13/19. 

□ 

Remark We have extended this method further, obtaining even smaller bounds for 
h, but as we increase the number of sites we consider, the calculations quickly become 
very elaborate, and the bound we obtain decreases very slowly. The smallest bound we 
recorded was 0.683. 

6. The growth rate for h < 0. 

In this section, we prove Theorem |3.4| The proof will follow closely the method used 
for the proof of Theorem 4.1 in 0, based on the estimates presented in [S], Section 2. 
An important ingredient to obtain the bounds is that after stabilization, every site has 
mass at most 1, so that the mass n starting from the origin, must have spread over a 
minimum number LfirsJ of sites. On the other hand, as shown in the proof of Theorem 
3.1, part 2, we know: \T\ < xrr- 

But crucial is the use of Green's functions and their asymptotic spherical symmetry, 
allowing to conclude more about the shape of the set of sites that split. Thus we can 
derive that T contains a ball of cardinality comparable to the coarse estimate [jzrhl ■ 
Moreover, T is contained in a ball of cardinality close to Li/frxJ- 

The method for ft, < does not depend on abelianness or monotonicity, therefore it 
can be adapted to the splitting model with arbitrary splitting order. 

We start with introducing some notation. 

Denote by Po,Eo as the probability and expectation operator corresponding to the 
Simple Random Walk {X{t)) starting from the origin. For d > 3, define 

oo 

g{z) = Eo^Ix{t)=z- 
t=i 

For d = 2, define 

n 

9n{z) = Eo^Ix(t)=z, 
t=l 

and 

g{z) = lim [g„(z) - g„(0)]. 

n— ^oo 

Defining the operator A as 

From [7], we have ^g{z) — — \ when x is the origin and ^g{z) = for all other x G U^. 

By u{x), we denote the total mass emitted from x during stabilization. Then Au(x) 
is the net increase of mass at site x during stabilization. For all x, let rjaoix) be the final 
mass at site x after stabilization. Since the final mass at each site is strictly less than 1, 
we have for all x G Z'' 

(6.1) lS.u(x) + (n — h)5Q^x — Voo{x) — h < 1 — h, 
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with Sq^x — 1 ii X is the origin and for all other x. 

Moreover, since the final mass at each site a; e T is in [0, 1), we have for all x € T 



(6.2) ^ h < Au{x) + (n — /i)(5o.a; = Vocix) — h < 1 ^ h, 

Proof for the inner bound: 

For X e Z'^, \x\ is the Euclidean distance from x to the origin. Let 
Ux) = (1 - h)\xf + {n- h)g{x) if d > 2, 

and let 

^d{x) = ld{x) - ld{\cir\ei), 

with ei = (1,0,0, ...,0). 

From Lemma 2.2 of [5], we have: 

U{x) = 0(1), X e 9Bei., 

therefore there is a constant C > such that |^d(a;)| < C, x G dUcir- Then 

u{x) - ^d{x) > -(,d{x) > -C, X e dBcir- 



Furthermore, using (6.1 1, we have for all x ^ Z'^, 



A(u - Cd) = Aw - ACd < 1 - - (n - h)So,x - (1 - /i) - (n - h)So,x = 0. 

Therefore, u — is superharmonic, which means that it reaches its minimum value on 
the boundary. Now, as in the proof of Theorem 4.1 of [8], the estimates in Lemmas 2.1 
and 2.3 of [S] can be applied to conclude that there is a suitable constant C2 such that 
u{x) is positive for all x € Bcir-C2- I— ' 

The proof for the outer bound is more involved that that in [8], because Lemma 4.2 
from [S], which is valid for the abelian sandpile growth model with ft, < 0, is not applicable 
for the splitting model. In essence, this lemma uses that if u{x) > u{y) for some x and y, 
then the difference must be at least 1, because mass travels in the form of integer grains. 
Clearly, we have no such lower bound in the splitting model. 

We note that in [5], a difi^erent proof for the outer bound appeared which is valid for 
the abelian sandpile growth model with h < d. Unfortunately, we cannot adapt this 
proof for the splitting model either. We will comment on this in Section [7[ 

We therefore first present some lemma's which we need to prove the outer bound. 

Lemma 6.1. Let h < 0, and take xq £ T adjacent to dT- 
There is a path xq xi ^ X2 ^ ■ ■ ■ ^ x,m — inT with 

u{xk+i) > u{xk) - ^^'^ ^ h, fc = 0, . . . , TO - 1. 



Proof. We will first show that we can find a nearest neighbor path such that: 

u(Xk+i) - u(Xk) + ^ 

Let Xi be the nearest neighbor of xq that loses the maximal amount of mass among all 
the nearest neighbors of Xq. If there is a tie, then we make an arbitrary choice. Because 
Xq has at least one neighbor that does not split, AM(a;o) < '^^^u{xi) —u{xo). Therefore: 
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(6.3) ~2d~^^^^^ ~ ~ ^""(^^o) = Vocixo) - h> -h. 

For fc > 1 take Xk+i ^ Xk~i to be the site that looses the maximal amount of mass 
among all the nearest neighbors of Xk except Xk-i- If there is a tie, then we make an 
arbitrary choice. It is always possible to choose Xk+i- As long as Xk is not the origin, 
then we get from (6.2) that: 

^'^^^(xfe+i) + ^^i^k-i) - u{xk) > r]oo{xk) - h> -h. 

Thus we get a chain {x^, fc = 0, 1, . . .} of nearest neighbors, that possibly ends at the 
origin. 

We rewrite 

2d , , 1 , , 2d , 



With this expression, and the fact that uq < ui by ( |6.3| , it is readily derived by induction 
that u{xk-i) < u{xk)- Inserting this, we obtain u{xk+i) > u{xk) — 2d-i ^^ ^'^ ^^^^ u{xk) 
is strictly increasing in fc (recall that /i < 0). 

Now it is left to show that the chain does end at the origin. We derive this by 
contradiction: suppose the chain does not visit the origin. Then the chain cannot end, 
because there is always a new nearest neighbor that loses the maximal amount of mass 
among all the new nearest neighbors. But the chain cannot revisit a site that is already 
in the chain, because u{xk) is strictly increasing in fc. But by (5.1), the chain cannot 
visit more than jj^^ sites. Therefore, the chain must visit the origin. □ 

Define Qk{x) = {y G Z"^ : max^ \xi — yi\ < k} a.s the cube centered at x with radius fc. 
Let 

iiW(x) = (2fc + l)-'^ ^ u{y) 
yeQk(x) 

be the average loss of mass of the sites in cube Qk{x), and 

r^"^ ^{x:Qk{x)cr}. 
Lemma 6.2. AuW(x) > - h - f^^loeQ.i.), for all x e T'^'^l 
Proof. From Proposition 5.3 of [6J, we know for every a:: 

(6.4) i]oo{y) > —(number of internal bounds in Qk{x)). 

yeQk{x) 



Equation (6.1 ) tells that Au{y) = rjady) — h — [n — h)6o^y. Therefore 

Au^^Hx) = [v^iy)-h-{n-h)So,y]. 
^ ' yeQk{x) 

Since Qk{x) has 2dk(2k + 1^^^ internal bounds, we get: 



2fc + 1 



(2fc + l)'^-^°^'5^(")- 



□ 
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Lemma 6.3. For every x ^ T*''^-', 

u{x) < a', 

where a' depends only on k, d and h. 

Proof. For x ^ T^^\ there is at least one site j/o G Qk{x) that does not split. For / > 1, 
take yi as the nearest neighbor of that loses the maximal amount of mass among all 
the neighbors of ?/;_!. Since j/g does not split, we have 



v~yo 

Therefore u{yi) < 2d{l — h). For every Z > 1, we have from (6.1) that ^ Tly^yi ^iv) < 
l — /i + u(y/), therefore < 2d{l — h) + 2du{yi). We know there are at most {2k + l)'^ 

sites in {?/;}/=o- Then: 

max u{x) <il-h) [(2d) + {2d)^ + ■■■ + (2d)(2'=+^'1 < 2(1 - h){2dY^''+^'^\ 

so we can choose a' = 2(1 - /^)(2c^)(2'=+l)^ □ 

Proof of the outer bound: 

First, we wish to find an upper bound for u{x) for all x wit h c\ r — 1 < < c'j^r, that 
does not depend on n. If x is not in T^''\ then we use Lemma 
For X e T^''\ take 



6.3 



For a fixed small e, we choose k such that 



k 1 

> - - 



2fc + 1 - 2 
For the fixed chosen fc, define 

M^)^ (2k + i)d My). 

Take 

Mx) = <t>d{x) - 4>d{lc'ir\ei). 

By calculation, we obtain A(j)d{x) = A(f)d{x) — 1/2 — e — /i — (n — ^)loeQfc(a;)- Then from 
Lemma 6.2 we know 

(6.5) A(w('=) - (j)d) = Au^^) - A<j)d > 0,Vx e T'-'^l 

This shows that u^'^-' — (/i^j is subharmonic on T*^*^'. So, it takes its maximal value on the 
boundary. We combine this information with some lemma's: 

• Lemma 2.4 of [S] gives that for all x, (j)d{x) > —a for some constant a depending 
only on d. 

gives that for every x € dT'^^\ u{x) < a'. 



6.3 



Lemma 

Finally, from Lemma 2.2 of jB], there is a £2 which only depends on e,d and h, 
such that for x with c[r — 1 < |a;| < c'^r, (j)d{x) < 02- 
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The first two lemma's imply that for x G dT^''\ u'^''\x) — (t)d{x) < a' + a, an upper 
bound that does not depend on n. Therefore, since u'-'^) — 0^ is subharmonic on T^''\ 

u{x) - (f)d{x) < a' + a,yx e . 

Combining this with the third lemma, we get: 

u{x) <C2 + a' + a,V.T G B^'^^ n T*". 

Therefore, there is an upper bound for u{x) that does not depend on n, for all x g 
Bc'^r ri T*' ■ From Lemma 6.3, we know also for x ^ , u{x) < a' . Summarizing all, we 
obtain that for all x with c'j^r — 1 < |a;| < c[r, u{x) < C, with C a constant that does not 
depend on n. 

To summarize, for all x with c[r — 1 < \x\ < c[r, u{x) < C, with C a constant that 
does not depend on n. 

Now it remains to show that a site that splits, must lie at a bounded distance from 
'Bc[r- This follows from Lemma 6.1 From every site Xq that splits, there is a path along 
which u{x) increases by an amount of at least — 2d~i ^ every step, and this path continues 
until the origin. Then along the way, this path must cross the boundary of B^'^ri and 

there u{x) < C. Therefore, we can choose = ^^^^557^- '-' 



7. Open problems 
Based on numerical simulations, we present some tantalizing open problems. 
7.1. A critical hi 



In Theorem [3^2] we give two regimes for h for which we know that the splitting model 
is explosive resp. robust. In between, there is a large interval for h where we can prove 
neither. We conjecture however that the two behaviors are separated by a single critical 
value of h, and that this value does not even depend on the splitting order. In dimension 
2, our simulations indicate that this critical h is 2/3. 

Conjecture 7.1. 

(1) For the splitting model on Z'^ , there exists a he — hc{d) such that for all h < h^ 
the model is robust, and for all h > h^, the model is explosive. 

(2) hc{2) = 2/3. 

7.2. The robust regime. 

We have proved Theorem |3.4| for all h < 0. We hoped to extend this result to all 
h < 1/2, by adapting the proof used in Section 3.1 of [5] for the abelian sandpile growth 
model (ASGM). However, the first step of this proof uses the fact that ii„ is nondecreasing 
in n, where u„ is the total number of topplings that each site performs in stabilizing 77^. 
This follows from abelianness of the topplings. Since the splitting model is not abelian, 
we were not able to adapt this proof to work for our model. Nevertheless, we conjecture 



Conjecture 7.2. Theorem 3.4 holds for all h < 1/2. 
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Figure 9. The splitting automaton with h = 0.667 and n ~ 16: From 
left to right, t = 17, t = 2A, t = 39, t = 76, t = 103. In this case, a 
limiting shape may not exist. 



7.3. The explosive regime. 

We have only just started classifying the multitude of shapes of the splitting automa- 
ton, that one can observe by varying h. We are confident that our method is capable of 
generating many more limiting shape results. In some cases, we observe that varying n 
can make a difference, however, we expect the following to be true: 

Conjecture 7.3. For the splitting automaton on , for every h E [1 — , 1) there 
exists a uq such that for every n > no, the limiting shape is a polygon, and depends only 
on h and d. 

This conjecture is reminiscent of Theorem 1 on threshold growth in 5 , but the splitting 
automaton is not equivalent to a two-state cellular automaton. 

For smaller values of h, the behavior of the splitting automaton seems to be not nearly 
as orderly. In Figure |9j we show the behavior at /i = 0.667, where we conjecture the 
model to be explosive. The shape of T seems to alternate between square and rounded. 
We are not sure whether a limiting shape exists for this value of h. 
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